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Abstract 

t  l  :  <<■'  ' 

^We  consider^  a  finite  collection  of  polyhedra  whose  defining  linear 

'  j 

systems  differ  only  in  their  right  hand  sides.  Jeroslow  [5]  and  Blair  [4] 
specified  conditions  under  which  the  convex  hull  of  the  union  of  these 
polyhedra  is  defined  by  a  system  whose  left  hand  side  is  the  common  left  hand 
side  of  the  individual  systems,  and  whose  right  hand  side  is  a  convex 

-(  '■  i  ‘  ,  ( 

combination  of  the  individual  right  hand  sides.  give"  a  new  sufficient 

condition  for  this  property  to  hold,  which  is  often  easier  to  recognize.  In 

it  >'* 

particular,  -wer  showfythat  the  condition  is  satisfied  for  polyhedra  whose 

defining  systems  involve  the  node-arc  incidence  matrices  of  directed  graphs, 

S'"  y  ~  ./•*• 

with  certain  right  hand  sides.  ^-We  also  derive^ as  a  special  case^  the  compact 
linear  characterization  of  the  two  terminal  Steiner  tree  polytope  given  in 

s 

Ball,  Liu  and  Pul leyb lank  — O  -  !  f- 
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1 .  The  Result 


Consider  a  collection  of  nonempty  polyhedra  in  Rn  of  the  form 
Pj  :=  {  x  e  Rn  J  A*x  =  d1 ,  x  i  0  },  ieT, 
where  T  is  a  finite  set  with  J T j  =  t  i  2,  and  for  i  e  T,  A1  is  an  m  x  n  matrix 
and  d1  is  an  m^ -vector. 

It  is  known  [1]  (see  [2]  for  a  published  version)  that  clconv  u  (P^ieT), 
the  (closed)  convex  hull  of  the  union  of  the  polyhedra  P^,  is  the  set  of 
points  x  e  Rn  that  have  an  extension  (x,  x1,  ....  xl,  A)  e  Rn+t<'*t  satisfying 
the  conditions 

x  -  E(x*:i  e  T)  =  0 

(1)  A‘xl  -  d‘\  =  i  e  T 

E(A^ : i  e  T)  =  1 

x* ,  A^  i  0,  i  e  T 

This  implies  that  in  all  basic  solutions  to  the  system  (1),  A^  =  0  or  1, 
i  e  T. 

The  question  naturally  arises,  is  there  a  more  compact  representation  of 
the  convex  hull  in  the  case  when  A1  =  A  for  all  i  e  T,  i.e.  when  the  polyhedra 
are  of  the  form 

PA  :=  {  x  e  R"  |  Ax  =  d1,  x  i  0  },  ieT. 

In  particular,  let  Q  be  the  set  of  those  x  e  Rn  that  have  an  extension 

(x.  A)  e  Rn*1  satisfying 

Ax  -  E(diA  :i  e  T)  =  0 

(2)  £(A4 : i  e  T)  =  1 
x  *  0,  Ax  i  0,  ieT, 

and  let  C  :=  cl  conv  u  (P  : i  e  T). 

It  is  easy  to  see  that  C  £  Q,  since  (2)  can  be  obtained  from  (1)  (when 

A  =  A,  i  e  T)  by  left-multiplying  the  first  equation  with  A,  and  then  adding 


1 


51 


to  it  all  remaining  equations  except  for  the  last  one.  However,  C  2  Q  is  not 
true  in  general. 

Jeroslow  [5]  gave  a  sufficient  condition  for  C  =  Q  to  be  true.  Blair  [4] 
gave  two  weaker  sufficient  conditions,  one  of  which  is  also  necessary  when  a 
certain  requirement  is  satisfied.  He  also  showed  that  recognizing  whether 
C  =  Q  is  NP-hard.  We  give  a  new  sufficient  condition  for  C  =  Q  which  is 
sometimes  easier  to  recognize. 

We  will  denote  by  d(X)  the  convex  combination  of  the  right  hand  sides  d* 
with  weights  X.,  i  e  T;  i.e.,  for  Xa  i  0,  i  =  1,  —  ,t,  such  that  ECX^i  e  T)  = 
1,  d(X)  =  E(diXi:  i  e  T) .  Also,  we  will  assume  that  A  is  of  full  row  rank. 

Theorem  1 .  C  =  0  if  for  every  a  x  m  noasingular  submatrix  B  of  A  and 
every  convex  combination  d(  X)  of  the  vectors  d  1 ,  i  e  T, 

(3)  B  1  d(  \)  i  0  implies  B  *d  1  i  0  for  all  i  e  T  such  that  X  >  0. 

i 

Proof.  Let  (1’)  be  the  system  obtained  from  (1)  by  replacing  each  A*  by 
A.  Since  C  £  Q  is  always  true,  we  have  to  show  only  that  Q  £  C  if  (3)  holds 
for  every  B  and  X.  Suppose  (3)  holds  and  let  x  e  Q.  W.-t.o.g.,  assume  (x,  X) 
is  a  basic  feasible  solution  to  (2).  Then  there  exists  an  m  x  m  nonsingular 
submatrix  B  of  A  such  that  x  =  (x  ,  0),  where  x  =  B  1d(X)  3t  0.  From  (3), 

0  B 

X,  >  0  implies  B  1d1  Jt  0.  Letting  x1  =  (x*,  0),  where  x*  =  B  *dAX  ,  we  have 

i  B  B  i 

that  (x,  x1,  ...,  x1,  X)  satisfies  (l’)»  hence  x  e  C.  || 

The  assumption  that  A  is  of  full  row  rank  is  not  essential.  If  A  does 
not  have  this  property,  we  replace  it  by  any  of  its  maximal  n-column 
submatrices  of  full  row  rank. 

2.  An  Application:  Unions  of  (Multiple)  Network  Polyhedra 

Next  we  discuss  an  important  class  of  problems  for  which  condition  (3)  of 
Theorem  1  is  always  satisfied. 

We  will  say  that  P  ,  i  e  T,  are  network  polyhedra.  if  A  is  the  node-arc 


incidence  natrix  of  a  directed  graph,  and  each  d1  has  one  component  equal  to 
some  positive  number  v,  one  component  equal  to  -v,  and  all  other  components 
equal  to  0.  If  v  =  1,  the  network  polyhedra  are  path  polyhedra.  The  P  , 
i  e  T,  will  be  called  multiple  network  polyhedra  (multiple  path  polyhedra)  if 
A  is  a  block-diagonal  matrix  with  each  diagonal  block  the  node-arc  incidence 
matrix  of  a  directed  graph,  and  each  d1  has  a  subvector  for  each  block  of  A, 
with  exactly  one  positive  and  one  negative  component  equal  in  absolute  value 
(equal  to  1  in  absolute  value),  all  remaining  components  being  0.  In  other 
words,  A  and  d1 ,  i  £  T,  are  of  the  form 


where  the  blanks  are  zero  matrices. 

Theorem  2.  C  =  Q  if  the  polyhedra  P  ,  i  e  T,  are  (multiple)  network 

polyhedra,  and  for  each  k  €  {1, ....  q} ,  either  the  negative  or  the  positive 

i  k 

component  of  d  is  in  the  same  position  for  all  i  €  T. 

The  proof  of  Theorem  2  will  use  the  following  auxiliary  result.  Let  A 

be  the  node-arc  incidence  matrix  of  a  directed  graph  Q  =  (V,E)  and  for  v  e  V, 
let  A^  be  the  matrix  obtained  from  A  by  deleting  row  v.  It  is  well  known  that 
if  G  is  (weakly)  connected,  the  rank  of  A  and  of  Ay  is  |v|  -  1,  and  that  any 
( | V | — 1 )  x  ( |  V | —1 )  nonsingular  submatrix  of  A^  can  be  made  lower  triangular  by 
appropriate  row  and  column  permutations.  In  addition,  every  such  matrix  has 
the  following  property. 

Lena  3.  Let  B  =  f  b ^  )  be  a  (  J  V\  -1)  x  (  |  ?\  -1)  lower  triangular 

submatrix  of  .  Then  B  1  =■  (  P  t  jJ  is  lower  triangular  and  for 

i  =  1,  ... ,  \  V\  -1,  j  <  i,  every  nonzero  0 .  has  the  same  sign  as  b^  ^  . 

Proof.  Let  the  rows  and  columns  of  B  be  indexed  by 

M  =  {l,...,m  =  |  V  |  —1} ,  and  for  R,  C  e  M,  let  be  the  submatrix  of  B  with 
rows  and  columns  in  R  and  C,  respectively. 


V»>  ■  '» 


,Wi 


Jh-Ai 


We  use  induction  on  a.  The  statement  is  trivially  true  for  a  =  1. 
Suppose  it  is  true  for  m  =  1,  —  ,k  and  let  m  =  k  +  1  i  2.  Then,  denoting 
K  =  {l,...,k  =  m-1} , 

B~1  _  Bk  i  0  __  Kl'l _ i_°_ 

nK  !  ,  ,  -1  nK,nK,-l  !  ,  -1 

B  ,  b  -b  B  (B  )  ,  b 

L  b  i  an  nn  n  K  bb 

By  the  induction  hypothesis,  for  i  =  1, _ ,k  every  nonzero  entry  in  row  i 

K  - 1 

of  (B  )  has  the  saae  sign  as  b  .  Hence  the  claim  of  the  Lemma  is  true  for 
K  i  i 

the  first  m-1  rows  of  B  On  the  other  hand,  the  i-th  entry  of  the  row 

K 

vector  B  ,  if  nonzero,  has  the  sign  opposite  to  that  of  b  (from  the  node-arc 

a  i  i 

incidence  relation).  Thus  -B  (B  )  £  0  and  therefore  every  nonzero  entry  in 

n  K 

row  m  of  B  1  has  the  same  sign  as  b'1,  hence  b  .  This  completes  the 

Bin  BB 

induction.  || 

k 

Proof  of  Theorem  2.  Since  each  matrix  A  ,  k  =  l,...,q,  has  a  redundant 

row,  we  delete  from  each  A  the  row  corresponding  to  the  unique  nonzero  entry 
i  k 

of  d  ■  which  is  in  the  same  position  for  each  i  e  T.  We  then  delete  the 

i  k 

corresponding  entry  from  each  d  ,  i  e  T;  and  if  for  some  k  e  {l,...,q}  the 

i  k 

remaining  nonzero  entry  of  each  d  ,  i  e  T,  is  negative,  we  also  change  the 

•  k  k  ^  k  | 

sign  of  A  and  of  each  d  ,  i  e  T.  Let  A  and  d  be  the  resulting  matrix  and 

vectors.  When  we  have  applied  this  operation  to  each  block  of  A,  the 

resulting  vectors  d  ,  i  c  T,  are  all  nonnegative  and  the  resulting  matrix  A  is 

of  full  row  rank.  The  polyhedra  P  are  now  defined  by  the  systems  Ax  =  d1 , 

x  £  0,  i  c  T. 

Let  m  be  the  number  of  rows  of  A.  Then  any  mxm  nonsingular  submatrix  B 
of  A  can  be  brought  by  row  and  column  permutations  to  the  form 


where  each  B^  in  a  square  lower  triangular  matrix  and  the  blanks  are  zeros. 


iMOOK. 
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Let  B  be  any  such  matrix;  then  its  inverse  is 


From  Lemma  3,  B~  is  lower  triangular  and  in  every  row  i  e  {l,...fm}  of 
B_1,  all  nonzero  entries  have  the  sign  of  b  .  Since  d1  i  0,  i  e  T,  for  any 

I  T  I  —  1 

X  e  R  ,  X  i  0,  such  that  ^(X^ii  e  T)  =  1,  B  d(X)  £  0  implies  that  every 
diagonal  element  bt  .  such  that  Xj  >  0  must  be  1  (i.e.  positive).  But  then 
B^d1  £  0  for  all  i  e  T  such  that  >  0,  i.e.  the  polyhedra  P.  satisfy 
condition  (3)  of  Theorem  1;  therefore  C  =  Q.|| 


3.  Another  Application:  Two  Terminal  Steiner  Tree  Polyhedra 
As  a  further  application,  we  dicuss  the  case  of  two  terminal  Steiner  tree 
polyhedra  and  derive  from  our  result  the  compact  linear  characterization 
obtained  (by  other  means)  by  Ball,  Liu  and  Pulleyblank  [3].  Let  G  =  (V,E)  be 

an  arc-weighted  directed  graph  with  three  distinguished  nodes:  a  source 

(root)  s,  and  two  sinks  (terminals),  t  and  r.  A  two  terminal  Steiner  tree  in 
G  is  a  minimum-weight  arborescence  rooted  at  s  and  containing  nodes  t  and  r. 
Such  an  arborescence  is  the  union  of  three  directed  paths,  from  s  to  some  node 
v,  from  v  to  t  and  from  v  to  r,  where  v  may  or  may  not  be  distinct  from  s,  t 
and  r.  However,  the  converse  is  not  true,  i.e.  the  union  of  directed  paths 
from  s  to  v,  from  v  to  t  and  from  v  to  r,  is  not  necessarily  a  rooted 

arborescence:  it  is  one  if  and  only  if  the  three  paths  are  node-disjoint. 

When  this  is  not  the  case,  the  union  of  the  three  paths  is  called  a  Steiner 


For  any  v  e  V,  the  incidence  vectors  of  directed  paths  from  s  to  v  in  G 
are  the  extreme  points  of  the  polyhedra  defined  by  the  system  Ax  =  e  -  e  , 

9  V 

x  i  0  where  A  is  the  node-arc  incidence  matrix  of  G  and  for  i  e  V,  e  is  the 

i 

I  V  I 

i-th  unit  vector  in  R  .  Also,  the  incidence  vectors  of  directed  paths  from 


*  ? 


ass* 


v  to  t  and  from  v  to  r  in  G  are  the  extreme  points  of  the  polyhedra  defined  by 
the  systems  Ax  =  e  -  e  ,  x  2t  0  and  Ax  =  e  -  e  .  x  i  0,  respectively. 

v  t  v  r 

I  E  I 

Now  let  P(v)  be  the  set  of  those  x  C  R  for  which  there  exist  vectors 
x*,  x*,  xr  e  R*E'  satisfying 

Ax’  =  e  -  e 

9  V 

(4)  Ax1  =  e  -  e 

v  t 


=  e  -  e 
v  r 


-lx"  -  lx1  -  Ixr  +  lx  =  0 


•  1  r  ^  n 

x  ,  x  ,  x  i  0. 


It  is  not  hard  to  see  that  the  extreme  points  of  P(v)  are  precisely  those 
x  e  R^E^  such  that  x  =  x’  +  xl  +  xr  for  some  incidence  vectors  x’,  xl  and  xr 
of  directed  paths  in  G  from  s  to  v,  from  v  to  t  and  from  v  to  r,  respectively. 
Indeed,  since  x  is  unconstrained  in  sign,  the  last  matrix  equation  of  (4)  does 
not  impose  any  constraint  on  x*,  x1,  and  xr .  Therefore,  x  is  an  extreme  point 
of  P(v)  if  and  only  if  each  of  x*.  xl ,  and  xr  is  an  extreme  point  of  the 
polyhedron  defined  by  x'  i  0  and  the  first,  xl  i  0  and  the  second,  xr  Je  0  and 
the  third  matrix  equation  of  (4),  respectively. 

To  get  from  (4)  a  system  of  full  row  rank,  we  delete  row  s  from  the  first 
matrix  equation,  row  t  from  the  second,  and  row  r  from  the  third  one,  and  we 
also  change  the  signs  in  the  first  matrix  equation.  We  thus  obtain  the  system 


-lx’  -  lx1  -  Ixr  +  lx  =  0 

*  1  r  ■*.  n 

X  ,  X  ,  X  i  0, 

equivalent  to  (4);  and  we  can  write  for  v  e  V, 


P(v)  = 


x  €  ^IeI  3  x*,  x1,  xr  c  R  E  such  that 
(x*,  x  ,  xr,  x)  satisfies  (5) 


i 

$ 

Wl* 


$ 
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Consider  now  the  two-terminal  Steiner  tree  problem  in  Q.  If  the  weights 
cj  are  nonnegative,  then  this  problem  is 

min  {  cx  |  x  £  conv  u  (P(v)  :  v  e  V)  }, 

since  the  minimum  is  always  attained  for  some  x  that  is  the  incidence  vector 

of  an  arborescence  rooted  at  s  and  continuing  nodes  t  and  r.  (If  the  costs 

are  arbitrary  and  there  exists  a  negative-cost  directed  cycle,  the  above 

problem  has  no  finite  minimum.  If  G  has  no  negative-cost  cycles  but  some  of 

the  costs  are  negative,  the  optimum  may  occur  for  a  Steiner  net  instead  of  a 

Steiner  tree,  and  some  arcs  may  have  a  flow  greater  than  one). 

Theorem  4.  The  ( closed )  convex  hull  of  the  union  of  P(  v)  for  all  v  e  V 

is  the  set  of  those  x  c  R^  for  which  there  exist  vectors  x  *,  a'  1 , 
r  I  E I  I V I 

x  e  R  and  X  e  R  satisfying  the  system 
-A  x’  -IX  =  0 

9 

(6)  A^  -IX  =  0 

A  xr  -IX  =  0 

r 

-lx*  -lx1  -  Ixr  +  lx  =0 

IX  =  1 

x*,  xl,  xr  Jt  0,  X  a  0. 

Proof.  The  system  (6)  can  be  rewritten  as 
My  -  LX  =  0 

(7)  IX  =  1 

y,  X  =  0 

where  L  =  (I,  I,  I)T  and  y  =  (x*,  x1,  xr ,  x)  e  R4IEL  We  will  show  that 

condition  (3)  of  Theorem  1  is  satisfied  for  (7)  viewed  as  an  instance  of  (2). 
The  matrix  M  has  3 1 V |  -  3  +  { B |  :=  p  rows  and  is  of  full  row  rank.  Any  p  x  p 
nonsingular  submatrix  of  M  can  be  brought  by  row  and  column  permutations  to 


the  form 


B  = 


J  J  J  I 

a  t  r 


where  B#,  B^  and  B^  are  (jvj  -  1)  x  (|v|  -  1)  lower  triangular  submatrices  of 
A#,  and  A^ ,  respectively,  J^,  and  are  { E (  x  ( | V J  -  1)  submatrices  of 
I,  I  itself  is  the  identity  matrix  of  order  jE|,  and  the  blanks  are  zero 


.-1 


matrices.  Clearly,  B  is  of  the  form 

- 1 


B 


-1 


B 


*  S'1 


1  B'1 

r 

*  *  * 
j  j  j  i 

sir 


A  A  A 

for  some  unspecified  J  ,  J  and  J  . 

at  r 

From  Lemma  3,  each  diagonal  block  of  B  1  is  lower  triangular  and  every 
nonzero  entry  in  any  of  the  first  3 | V | —3  rows  of  B_1  has  the  same  sign  as  the 
corresponding  diagonal  element  of  B.  It  then  follows  that  for  B_1LX  £  0  to  be 
satisfied  for  any  X  i  0  (with  IX  =  1),  it  is  necessary  that  every  diagonal 
entry  of  B  be  +1.  But  then  B  *LX  £  0  implies  B  'f-1  £  0  for  every  column  -£*  of 
L  such  that  X^  >  0,  i.e.  condition  (3)  of  Theorem  1  is  satisfied. || 

Theorem  4  asserts  that  all  basic  solutions  to  the  system  (6)  are  integer, 
and  thus  (6)  represents  a  compact  (polynomial-sized)  linear  characterization 
of  the  two  terminal  Steiner  tree  polytope.  This  characterization  is  not  new; 
in  fact,  it  is  due  to  Ball,  Liu  and  Pulleyblank[3] .  However,  our  proof  of  its 
validity  is  new,  and  puts  this  problem  in+o  the  more  general  context  of  unions 
of  polyhedra  whose  convex  hull  has  the  simple  representation  (2). 

It  should  be  mentioned  that  the  approach  discussed  here  can  be  extended 
to  k-terminal  Steiner  tree  polytopes;  the  resulting  formulation  involves  the 
union  of  a  number  of  polytopes  exponential  in  k,  but  polynomial  in  |v|  for 
fixed  k. 
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